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I. INTRODUCTION 



The vector (V) and axial-vector (A) current-current correlators are fundamental quantities of the strong-interaction 
physics, sensitive to small- and large-distance dynamics. They serve as an important testing ground for QCD as well 
as for effective models of strong interactions. In the limit of the exact isospin symmetry the V and A correlators in 
the momentum space (with — q 2 = Q 2 > 0) are defined as 

g^q 2 )Ti v AQ 2 )5 a \ (1.1) 
9a V q 2 ) Il$(Q 2 )S ab + q^Ut(Q 2 )6 ab , (1.2) 

where the QCD currents are 

J» = q^T a q, Jl a = q^ l5 T a q, (1.3) 

and T a denote the generators of the SUf(2) flavor group, normalized to txT a T = \8 ah . The momentum-space 
two-point correlation functions obey a (suitably subtracted) dispersion relation, 

n(Q2) = -J„ 7T^ Imnw - (L4) 

where the imaginary parts of the correlators determine the spectral functions 

ui(s) = 47rlmn^(s + iO), 

ffll (s) = 47rlmn^(s + i0). (1.5) 



nl'«\q) = i J d'x e^Ulfix) = (q,q„ - 

n#*(*) = <o|T{j>)j>)t}|o), 
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Recently, the inclusive nonstrange V and A spectral functions have been determined separately and with high precision 
by the ALEPH and OPAL collaborations from the hadronic r-lepton decays (r — > v T + hadrons) in the interval 
of invariant masses up to the r mass, < s < m 2 . 

The difference of the V and A correlation functions is very sensitive to the details of the spontaneous breaking of 
the chiral symmetry. In particular, the behavior of this combination is constrained by the chiral symmetry in the 
form of sum rules obtained through the use of the optical theorem 0. [H 0- . The experimental separation of the 
V and A spectral functions allows us to accurately test the chiral sum rules in the measured interval 

013 The 

coefficients of the Taylor expansion of the correlators at Q 2 — are expressed in terms of the low-energy constants of 
the chiral perturbation theory(xPT) Q . On the other hand, the large-s behavior of the correlators can be confronted 
with perturbative QCD thanks to the sufficiently large value of the r mass. In the high-s limit H v (Q 2 ) and H A (Q 2 ) 
arc dominated by the free-field correlator, corrected by nonperturbative terms with inverse powers of Q 2 . This follows 
from the fact that the correlators can be represented by an operator product expansion (OPE) series and thus are 
sensitive to the nonperturbative physics at smaller energy scales Q • Recently, the interest in the OPE expansion has 
been revived due to a possible appearance of unconventional quadratic power corrections, ~ 1/Q 2 , found in @, and 
also observed in lattice simulations Q. The ALEPH and OPAL data have been intensely used in the literature in 
order to place limits on the leading coefficients of the xPT and OPE expansions (see, e.g, [Tol ITU IT^. Il3l Hil ITEf ) . 

The aim of this work is to calculate the nonperturbative parts of the V and A current-current correlators in the 
kinematic region reaching up to moderately large Euclidean Q 2 and to extract experimentally observed characteristics. 
The calculations are carried out in the effective chiral model with nonlocal quark-quark interactions, which is made 
covariant by the inclusion of the P-exponents in the non-local interaction vertex. A specific prescription for the Wilson 
lines and their differentiation, described in Sec. Ill, follows exactly Ref. 0]. That way the model is made consistent 
with the gauge invariance and can be used to analyze the V and A correlators. The model is a nonlocal extension 
of the well-known Nambu-Jona Lasinio model. Moreover, its nonlocal structure may be motivated by fundamental 
QCD interactions induced by the instanton and gluon exchanges, which induce the spontaneous breaking of the chiral 
symmetry and generate dynamically a momentum-dependent quark mass. From the point of view of the standard 
OPE, the whole series of power corrections characterizes nonlocal properties of the QCD vacuum and may be described 
in terms of the nonlocal vacuum condensates |l7l ITsj . The use of a covariant nonlocal low-energy quark model based 
on the self-consistent approach to the dynamics of quarks has many attractive features, as it preserves the gauge 
invariance, is consistent with the low-energy theorems, as well as takes into account the large-distance dynamics 
controlled by the bound states. Similar models with nonlocal four-quark interaction have been considered earlier in, 

e.g., m m m m sin mm mm 

and applied to describe various low-energy phenomena. 

Nonlocal models have an important feature which makes them advantageous over the local models, such as the 
original Nambu-Jona-Lasinio model. At high virtualities the quark propagator and the vertex functions of the quark 
coupled to external fields reduce down to the free quark propagator and to local, point-like couplings. This property 
allows us to straightforwardly reproduce the lea ding terms of the operator product expansion. For instance, the second 
Weinberg sum rule is reproduced in the model [25(, which has not been the case of the local approaches. In addition, 
the intrinsic nonlocalities, inherent to the model, generate unconventional power and exponential corrections which 
have the same character as found in Q and in the instanton model (see e.g. |l3j|). Recently, the nonlocal effective 
model was successively applied to the description of the data from the CLEO collaboration on the pion transition 
form factor in the interval of the space-like momentum transfer squared up to 8 GeV 2 p9l l30| . Importantly, in that 
study at zero photon virtualities the chiral anomaly result were reproduced, while at high photon virtualities the 
factorization of short and long distances occurs at a scale of the order of 1 GeV 2 . This allowed for the extraction of 
the pion distribution amplitudes of leading and next-to-leading twists. There are several further advantages in using 
the nonlocal models compared to the loc al ap proaches. The non-local interactions regularize the theory in such a 
way that the anomalies are preserved [3ll l32j |. In other regularization methods in the local models [331 I'll l35| the 
preservation of the anomalies can only be achieved if the (finite) anomalous part of the action is left unregularized, 
and only the non-anomalous (infinite) part is regularized. Next, with non-local interactions the model is finite to 
all orders in the 1/N C (loop) expansion. Finally, as shown in |27|. stable solitons exist in a chiral quark model with 
non-local interactions without the extra constraint that forces the chiral fields to lie on the chiral circle. 

In the present paper we further test the nonlocal quark model by carrying out an analysis of the momentum de- 
pendence of the current-current correlators. A transformation of the spectral functions measured by the ALEPH 
collaboration into the Euclidean momentum space allows us for a precise an unambiguous comparison of the experi- 
mental data with the model calculations. The paper is organized as follows. In Sect. II, we briefly recall the results 
of the chiral perturbation theory and operator product expansion concerning the V and A correlators. In Sect. Ill 
and IV, we outline the gauged nonlocal quark model and introduce the quark-current vertices. Then, we derive the 
expressions for the nonperturbative parts of transverse V and A correlators (Sect. V) and, after fixing the model 
parameters in Sect. VI, confront the results with the available experimental data at large (Sect. VII) and low (Sect. 
VIII) Euclidean momenta. We explicitly demonstrate the transverse character of the V and nonsinglet A correlators 
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in Sect. IX. In Sect. X the contribution of the ?7a(1) axial anomaly to the flavor-singlet longitudinal A correlator 
is displayed and the topological susceptibility is calculated as a function of the momentum. Below, in all cases we 
use the strict chiral limit, with current quark mass equal to zero. Numerically current quark mass corrections to 
observables discussed in the paper do not exceed more than few percents. 

II. CHIRAL SUM RULES AND THE OPERATOR PRODUCT EXPANSION 

Chiral sum rules are dispersion relations between the real and absorptive parts of a two-point correlation function 
that transforms symmetrically under SU(2) l x SU(2)r (for the case of non-strange currents). Through the use of the 
dispersion relations the sum rules are directly expressed in terms of the difference of the V and A spectral densities. 
Here is the list of sum rules, given in the strict chiral limit, which are investigated in this paper: 

The first Weinberg sum rule (WSR I) 0, 

4^2 l ds [ Vl (s) - fll («)] = [-Q 2 I%~ A (Q 2 )] Q2 ^ = fl (2.1) 
the second Weinberg sum rule (WSR II) 0, 

1 i-so-^oo 

j-S / dss [in (s) - a, (s)} = Q 2 [-Q 2 T%- A (Q 2 )] Q2 ^ = 0, (2.2) 

•J 

the Das-Mathur-Okubo (DMO) sum rule 0, 

1 rsa^oo i q 

I [vx (s) ~ ax «] = ^ [Q 2 U V T - A (Q 2 )] 

and, finally, the Das-Guralnik-Mathur-Low-Yuong (DGMLY) sum rule p|, 

1 rsa^roo roc Aqr f 2 

-^J o dss In - [vi (s) - at (s)] = J dQ 2 [~Q 2 U^ A (Q 2 )} = [m 2 ± - m 2 ] , (2.4) 

where in the last equation a w 1/137 is the fine structure constant. In the chiral limit of massless quarks the DGMLY 
sum rule is independent of the arbitrary normalization scale, ji 2 , thanks to WSR II. It was shown by Witten 36] that 
the positive electromagnetic mass shift of the charged pions is a consequence of the DGMLY sum rule combined with 
the positivity property of the V — A correlator, 

- Q 2 U^,- A (Q 2 ) Ss 0, for < Q 2 < oo. (2.5) 

According to Witten, if the bare u and d quarks were massless and the mass shift were negative, the charged pions 
would become tachyons destabilizing the QCD vacuum. 

Whereas WSR I and DMO are low-energy sum rules (in the sense that the right-hand sides involve correlators at 
low momenta), and are reproduced in most low-energy effective quark models, WSR II is a high- momentum sum rule. 
In local models it is not reproduced, as discussed shortly. The DMGLY sum rule collects contributions from the whole 
range of Q 2 , both soft and hard. 

The left hand sides of the sum rules Ij2.1|l - H2.4(l have been determined with the experimental data from Q and 0, 
with so taken as the upper limit of the interval of the invariant mass covered by the experiment. The right-hand sides 
of the sum rules are the theoretical predictions as so — > oo. The DMO sum rule relates the derivative of Q 2 times the 
correlator to the square of the pion decay constant f n = (92 A ± 0.3) MeV [3]} obtained from the decays it~ — » iiTv^ 
and tt" — > /i - 77 M 7, to the pion axial-vector form factor Fa = 0.0058 ± 0.0008 for the radiative decays tt~ —* i _ P/7, 
and to the pion charge radius-squared (r 2 ) = (0.439 ± 0.008) fm 2 obtained from a one-parameter fit to the space-like 
data |38|. 

The listed chiral sum rules provide important restrictions on the correlators at low and high energies. The first 
Weinberg sum rule (|2.1|) fixes the normalization of correlators and holds in all variants of the Nambu-Jona-Lasinio 
models, local or nonlocal. In general, the coefficients of the Taylor expansion of the correlators at low Euclidean 
momenta are given by the low energy constants of the strong chiral Lagrangian. The second Weinberg sum rule H2.2(l 
signals that the leading asymptotics of the high Q 2 power expansion of the V — A correlator essentially starts from 
dimension d = 6 term, and as such is valid in the nonlocal versions of the effective chiral quark models |25J . In local 



= f 



F A 



(2.3) 
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models WSR II involves on the right-hand side the large constituent quark mass times quark condensate, M q < qq >, 
thus is violated badly. In this regard the nonlocal models are highly rewarding. 

More detailed short-distance, or large Q 2 , properties of the correlators are represented by the QCD operator product 
expansion @|. For the V — A and V + A combinations the OPE provides the following leading-twist terms in the 



chiral limit: 



n^(Q 2 ) - E % i = % £ + ^)' ( 2 - 6 ) 



d=6,8.. 



h^(q 2 ) = J2 = ( 2 - 7 ) 



d=0,2,4.. 



4tt 2 V W% 2 4tt 3 Q 2 12 Q 4 Q 6 l Q 8 



where the vacuum matrix elements of dimension d = 6 operators are 

O v „ A = na s [((u^X a d) (d^X a u)) - ((u W X a d) (dj^ 5 X a u))] , (2.8) 



Oy +A = -ira s 



((u^ 5 \ a d) (d W \ a u)) + {(u^\ a d) (d^X a u)) + \ J2 E ((Wilfj.X a qi) ijljl p,X a qj)) 



i—u^dj—u^d^, 



with A a being the color SU (3) matrices. 

The V — A correlator does not acquire any perturbative contribution in the limit of massless quarks, hence it 
is sensitive entirely to the chiral symmetry breaking parameters. Already at relatively small Q 2 the d = 6 term 
dominates in the expansion of IIy_^. On other hand, the sum of the correlators, TLy+A, supplied with small power 
corrections, is close to the free-field result for distances up to 1 fm In the expansion of HV +J 4 we also included 
d = 2 term which violates the original OPE expansion of 0- Motivation for inclusion this term into consideration 
was given in In (|2.6(l and (|2.7(l we do not explicitly show exponentially suppressed terms that may be induced 
by existence of instantons [?J . The magnitudes of the vacuum matrix elements which appear in the OPE can not be 
fixed from the first principles and are fitted to various hadronic observables. Clearly, their determination is bound to 
carry experimental and theoretical uncertainties. Only a sign of the d = 6 term in the V — A correlator is fixed by 
the Witten inequality: O v _ A < 0. 

Different models are used to estimate the vacuum expectation values. The standard approach in the calculation 
of the dimension d = 6 matrix elements suggested and used in original work was to explore the factorization 
hypothesis, i.e. the saturation of the four-quark matrix elements with the intermediate vacuum state. Under this 
assumptions for the dimension d = 6 matrix elements one gets 

[0 6 V . A ] faCt ° r = -32vra s (qqf /9, [O v+A ] ^ = 64™, (qq? /81. (2.9) 

However, some authors conclude that probably the factorization hypothesis is violated by a factor of 2-3 [33. Moreover, 
quite different result appears if one uses the instanton liquid model to calculate these matrix elements |40l l4l) . 

In this work, for comparison with other model results, we use the following typical values of the condensates found 
via standard QCD sum rules without and with the inclusion of the d = 2 term: 

— A 2 = 0, (—(G a u S) =0.012 GeV 4 , a s (qq) 2 = 2.4 • 1(T 4 GeV 6 [7,11], (2.10) 

^■A 2 = —0.12 GeV 2 , ^— (G°„) 2 ^> = 0.022 GeV 4 , a s (qq) 2 = 5.8 • 10~ 4 GeV 6 [42]. (2.11) 

The above values hold at a typical renormalization scale of about 1 GeV. 



III. GAUGING NONLOCAL MODELS 



In local theories, the gauge principle of the minimum action uniquely determines the interaction of the matter fields 
with the gauge fields. However, in nonlocal theories such an interaction may be introduced in various ways, and its 
transverse part cannot be uniquely defined |43| . In order to obtain the nonlocal action in a gauge- invariant form with 



FIG. 1: Diagrammatic representation of the effective nonlocal four-quark interaction of Eq. (13.21 . The hatched blobs represent 
the nonlocal interactions, with the Vi matrices present, and the dotted line indicates that the diagram can be cut across this 
line without crossing the quark lines going across the blobs. 



respect to external fields V and A, wc define the delocalized quark field, Q, with the help of the Schwinger gauge 
phase factor, a.k.a. the Wilson line or the link operator, 

Q(z, y) = Pexp U jf" dz^ [V%{z) + A<^(z) l5 ] T a J q(y), Q(x, y) = Q\x, y)^ . (3.1) 

Here V£(z) and A^(z) are the external gauge vector and axial- vector fields, respectively, and P is the operator of 
ordering along the integration path, with y denoting the position of the quark and x being an arbitrary reference 
point. The P operator arranges the matrices in each term of the expansion of the exponent from the left to the right 
in the order determined by the point z moving along the path from x to y. 

We start with the nonlocal chirally invariant action which describes the interaction of soft quark fields. The nonlocal 
four-quark interaction is depicted in Fig. ^ The soft gluon fields have been integrated out. The corresponding gauge- 
invariant action for quarks interacting through nonperturbative exchanges can be expressed in a form similar to the 
Nambu-Jona-Lasinio model 16] 



d^x 5(1)7" [id^ - y M (x) - 75 A M (a?)] q(x) + 
+ \g J d 4 X J Y[^x n f(x n ) \Q(X - Xl ,X)TiQ(X,X + x 3 )Q{X - x 2 ,X)riQ{X,X + x 4 )] , (3.2) 

where in the simplest version of the model the spin-flavor structure of the interaction is given by matrix product 

(r^Ti) = (1® l-H75T a ®H5T a ). (3.3) 

In Eq. (|3.2|) q = (u, d) denotes the quark flavor doublet field, G is the four-quark coupling constant, and r a are the 
Pauli isospin matrices. 

The derealization of the quark fields with the inclusion of the path-ordered Schwinger phase factors, Eq. (|3.1fl . 
ensures the gauge invariance of the nonlocal action 13.211 . However, the presence of these factors modifies the quark- 
current interaction, as shown graphically in Fig. J2J. The modification of interaction, required by the gauge principle, 
poses a technical difficulty in dealing with nonlocal models, as many diagrams appear in the analysis of physical 
processes. The ambiguities in making the nonlocal 4-quark interaction gauge invariant are manifest in the path- 
dependence in the definition Ij3.1|l . as well as in the choice of the junction of the quark sources with the gauge strings. 
In general, the Noether currents consist of two components: the path- independent longitudinal part and the path- 
dependent transverse part. The dependence of the transverse component on the choice of the path is a feature of 
any method employed in constructing the Noether currents corresponding to a nonlocal action, and this freedom is 
immanent to the formulation of the model. Wc should recall here that the discussed ambiguities in the construction 
of the transverse parts of the Noether currents are by no means specific to chiral quark models. They also appear, 
e.g., in nuclear physics when one considers meson-exchange processes. To summarize, the choice of the path in Eq. 
(|3.1|) is a part of the model building. 

In what follows, we use the formalism [23, ^| based on the path-independent definition of the derivative of the 
integral over a line for an arbitrary function F^z): 



J* dz v F v {z) = F^y), 6^ (x - y) £ dz v F v {z) = 0. 



(3.4) 



This choice effectively means that the differentiation involves moving the end-point of the line only, with the other 
part of the line kept fixed. As a result, the terms with nonminimal coupling, which are induced by the kinetic term 
of the action, are omitted. 
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FIG. 2: Diagrammatic representation of gauging of the effective nonlocal four-quark interaction 13.21 . shown up to the second 
order in the expansion in external fields, represented by the wavy lines. Note that vertices involving one quark line and multiple 
gauge fields appear. 

In general, external fields entering into Eq. (|3.2() are arbitrary auxiliary fields; however, some of them can be 
associated with electromagnetic, weak, or strong interactions. In the case of the electromagnetic interactions, the 
gauge factor takes into account the effects of the radiation of the photon field when the two quarks are moving 
apart. This formalism was used in pol I2H I23L I24I l25j to represent the nonlocal interaction in a gauge- invariant 
form. The incorporation of a gauge-invariant interaction with gauge fields is of principal importance if one desires to 
treat correctly the hadron characteristics probed by external currents, such as hadron electromagnetic and weak form 
factors, structure functions, distribution amplitudes, etc. 

In Eq. I|3.2|) the functions f(x n ), normalized to /(0) = 1, form the kernel of the four-quark interaction and 
characterize the space-dependence of the order parameter of the spontaneous chiral-symmetry breaking. Thus, the 
interaction is treated in the separable approximation. The choice of the nonlocal kernel in the form of (|3.2[1 is motivated 
by the instanton- induced nonlocal quark-quark interaction 0, where the nonlocal function f(x n ) is related to the 
quark zero mode emerging in the instanton field [TH Il9j |. To have the same flavor symmetry as in the original 
instanton-induced 't Hooft determinant interaction one needs to add yet another piece of the form 

G' (r a (g>T a + z 75 ® i75 ) , (3.5) 

with the coupling G' = — G. This term will be important in the discussion of the isosinglet axial currents in Sect. X. In 
the present work we do not consider an extended version of the model that explicitly includes vector and axial- vector 
degrees of freedom [23| (we take Gy = 0, therefore gA = 1 and M v , M\ — > 00). 

In order to compute physical quantities we must first determine the full quark propagator and the full vertices 
for the vector and axial- vector currents. All calculations will be done in the leading order of the 1/N C expansion, 
also referred to as the one-quark-loop level or the ladder approximation. In the nonlocal model the dressed quark 
propagator, S(p), with the momentum-dependent quark scalar self-energy (mass), M(p), is defined as 

S- 1 (p)=p-M(p). (3.6) 

Note that the considered model involves a constant quark wave-function renormalization function, Z(p) = 1. The 
equation for the quark propagator in the ladder approximation, also known as the gap equation, 

M M - UGN,N,f( P ) J ^f {k)w M^_ (3.7) 
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has the formal solution [2j| of the form 

M(p) = M q f(p), (3.8) 

with constant M q = M(0) determined dynamically from Eq. I|3.7jl . The quark self-energy is depicted in Fig. 10). Note 
that the functions f(p) are treated non-dynamically, i.e. their dependence on p is fixed, while M(p) is dynamical. 
Furthermore, the integrals over the momentum are calculated by transforming the integration variables into the 
Euclidean space, (k° —* ik^, k 2 — > —k 2 ). 



IV. CONSERVED VECTOR AND AXIAL- VECTOR CURRENTS 



The Noether currents and the corresponding vertices are formally obtained as functional derivatives of the action 
(|3.2|) with respect to the external fields at the zero value of the fields. For our purpose, it is necessary to construct the 
quark-current vertices that involve one or two currents (contact terms) . In the presence of the nonlocal interaction 
the conserved currents include both local and nonlocal terms. In order to expand the path-ordered exponent in the 
external fields, we use the technique described in [2(j (see also 0,H^|). Briefly, this method is as follows. First, the 
Fourier transform of the interaction kernel in Eq. (|3.2(1 is obtained and expanded in the Taylor series in momemta. 
Next, the momentum powers are replaced by the derivatives acting on both the path-ordered exponent and the quark 
fields. Finally, the inverse Fourier transform is performed and summation is carried out again. The relations H3.4(l 
and 



d 4 xF (x 2 ) e~ ipx I + dz" e- iqz = iX (2p + qXf F (P + X<1 \ F ^ e~ lqa + F (p) [ dz^ e~ lqz , (4.1) 
J y (p + M) -P 2 Jy 



where F(z 2 ) is an arbitrary function, are frequently used in the procedure described above |72j . The longitudinal 
projection of the above relation is 



d i xF (x 2 ) e- ipx J dz 11 =i[F(p+ Xq) e -' lqa - F (p) e~ i?y ] . (4.2) 

The algebra needed to obtain the vertices with this method is straightforward but somewhat tedious, hence below we 
present only the final results. 

The vector vertex following from the model ll-i.2t is (Fig. 

Tl(k, <?, k' = k + q)= T a [ 7|t - (k + fc')] , (4-3) 

where M^(k, k') is the finite-difference derivative of the dynamical quark mass, q is the momentum corresponding to 
the current, and k (k') is the incoming (outgoing) momentum of the quark, k' = k + q. The finite-difference derivative 
of an arbitrary function F is defined as 

(4 .4, 




FIG. 3: Diagrammatic representation of the quark scalar self-energy of Eq. 13.71 . 




Thus, with the gauging prescription given by (|3.2[l and l|3.4|) . one gets the minimum-coupling vector vertex without 
extra transverse pieces. The form of the vertex is the same as the longitudinal vector vertex resulting from the 
Pagels-Stokar construction [48|. The vertex satisfies the proper Ward-Takahashi identity: 



q^Jk, g, k') = Sp 1 (k') T a - T-Sp 1 (k) . 



(4.5) 



The vector vertex (14 .3|) is free of kinematic singularities and for this reason was advocated long ago in [43, |43 ■ For 
the case of the momentum-independent mass, as in local models, the vertex (|4.3|) reduces to the usual local form, 

The bare axial-vector vertex obtained from the action l|3.2[) by the differentiation with respect to the fields is given 
by the formula (c/. Fig. ^} 



Tl a {k,q,k' = k + q) = T a 



1, - (k + k%± 



k' 2 - k 2 



(4.6) 



-%2 v / M(fc')M(fc) 



75, 



where we have introduced the notation 

J A p(q 2 ) 

JAp{q 2 



M 2 

0) = -^-q 2 J' AP (0) + O(Q i ), 



with 



J' AP (0) = 



N c Nf f uM (u) [AM' (u) + 2uM" (u)] - u (M' (u)) 2 



32tt 2 



D(u) 



where u — k 2 and (in the Euclidean space) 



D(k) = /c 2 + M 2 (fc). 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



In Refs. [la, |45j it was demonstrated that in order to obtain the full vertex corresponding to the conserved axial- 
vector current it is necessary to add the term which contains the pion propagator. The presence of this term is 
associated with the well-known pion-axial vector mixing. The addition of the term with the pion propagator exactly 
cancels the third term in Eq. I|4.6[) . and the full conserved vertex acquires the form (c/. Figs. © and ©) 



r*?(k,q,« = k + q)=T a 



k ,2_ k 2 (y/MWj ~ VM(fe) 
1 ~ )m — 



k' 2 - k 2 



75- 



(4.11) 



It satisfies the axial Ward-Takahashi identity, 

q^ik, q, k') = j.Sp 1 (fc+) T a + T-Sp 1 (fc_) 75 . 



(4.12) 
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The axial- vector vertex has a kinematic pole at q 2 = 0, a property that follows from the spontaneous breaking of the 
chiral symmetry in the limit of massless u and d quarks. Evidently, this pole corresponds to the massless Goldstone 
pion. 

We also need the vertices that couple two currents to the quark (cf. Fig. In this regard it is convenient to 
introduce the notation 



Gl (k,q)= i T a (2k + q) u f( 1 \k,k + q), 



and 



Gf v (k,q,q',k') = -f(k'){T a T b [g^f^(k 1 k + q + q , )+ 

+ [2 (k + q') + g] p (2k + q') v (k,k + q',k + q + q')] + [(q, a, M ) «— (q' , b, v)}} , 
where the second finite-difference derivative is defined by 

Further, we need to introduce 

F^ a (k,q) = G«(k,q)±G*(k-q,q), 
F^ ab (k,q,q',k') = Gf v {k,q,q',k') + G ab v {k' - q - q' ,q,q' ,k) 



(4.13) 
(4.14) 



± G a (k, q) Gl (k' - q', q') ± G a (k' - q, q) G h v (k, q') . 



With the above definitions one gets for the VV contact term 

Tf v {k,q,q',k' = k + q + q') = M v F™ ab (k,q,q> \k' = k + q + q') 



f(k)f(k')G 



dH 

(27T) 4 



rTr 



S(l)F£ )ab (l,q,q',l) 



For the AA contact term one finds 



where 



T% ab (k, q,q',k' = k + q + q 1 ) = T 5 { ^(k, q,q',k' = k + q + q') + Ar^ b (fc, q,q',k l = k + q + q') 

T^l u (k,q, q \k' = k + q + q') = M q F^ ab (k, q, q' , k' = k + q + q') 

d 4 l 



f(k)f(k')G 



(2n) 4 



S(l)F^ ab (l >q ,q'J) 



An additional contribution appears for the AA iso-triplet contact term 

ATlf(k, q,q',k> = k + q + q') = —G [r c G* (fc, q') -G b u (k~ q', q') r 

dH 



(2-) 4 



rTr 



S(l)r c F^ a (l,q) 



+ [(9,M,a) < — > (q',v,b)] ■ 



(4.15) 

(4.16) 
(4.17) 

(4.18) 



(4.19) 

(4.20) 
(4.21) 

(4.22) 




FIG. 5: Diagrammatic representation of the full axial-vector vertex obtained from the bare axial vertex supplied with rescat- 
tering process of Fig. HJ In the present model there is no rescattering in the vector channel, and Fy = IV. 
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FIG. 6: Diagrammatic representation of the quark-quark scattering matrix, T. 



In the above expressions Tr denotes the trace over flavor, color, and Dirac indices. 

In the following we also need to introduce the polarization operator in the pseudoscalar channel (cf. Fig. 0), 

Jpp(q 2 )S ab = -^ f -^M(k)M(k + q )Tr[S(k) l5 T a S(k + q) l5 T b ] 
lvl q J (2n) 

and the correlator of the axial current vertex l|4.6fl and the pion vertex (cf. Fig. |SJ) 

T£ (k,k') =^75/(fc)/(fc')r Q , 

defined by 



(4.23) 



(4.24) 



J ^ {q2)Sab=9 fJ-0f Tr 



S(k)T^(k, q, k + q)S (k + q) T° (k + q, k) 



(4.25) 



Through the use of the gap equation (|3.7[) and the expression for the pion decay constant, / w , given by |l9l 12 



M(uf - uM(u)M'(u) + u 2 M'(uf 



D 2 (u) 



(4.26) 



these correlators have the following expansion at zero momentum: 



■Mi 2 ) = \ + %f + ° ^ ' J * A ^ = f2 + ■ 



(4.27) 



In i|4.26|l we have used the notation u = k 2 and M'(u) = dM(u) / du. In Eq. (|4.24|) the quark-pion coupling, 
gl = [J' PP (0))~ , and the pion decay constant, /„■, are connected by the Goldberger-Treiman relation, 



(4.28) 



which is verified to be valid in the nonlocal model [23l | , as required by the chiral symmetry. 



V. CURRENT-CURRENT CORRELATORS (TRANSVERSE PARTS) 



Our goal is to obtain the nonperturbative parts of the current-current correlators from the effective model and to 
compare them with the existing r decay data. The current-current correlators may be represented as a sum of two 
terms, dispersive (Fig. |5J and contact (Fig. I10[l . namely 



Q 2 K (Q 2 ) = K'l (Q 2 ) + (Q 2 ) 



(5.1) 




FIG. 7: The polarization operator of Eq. 14.231 
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FIG. 8: The correlator of the bare axial current vertex 114.611 and the pion vertex I4.24H . 



K> (Q 2 ) = / 7TT4 Tr K ( fc < Q,k + Q)S(k + Q) T v v (k + Q, -Q, k) S (*)] 

J (27T) 



K (Q 2 



(27T) 
(i 4 fc 

(2*)< 



-Tr 



(fc, Q,fc + Q)5(/c + Q) (fc + Q, -Q, k) S {k) 



S J ^(Q 2 ) = 2M q j A±.Tr[S(k)T J ^(KQ,-Q,k + Q)} 



(5.2) 
(5.3) 
(5.4) 



The vertices (k, q, k') are given in Eqs. (|4.3l4.11j) . T^ v (k, q, q' , k') in Eq. (j4.18l4.19jl . and T* u (fc, 9, g', fc') in Eq. 
Q4.6J1 . The difference in the definitions of (Q 2 ) in (|5.2(l and (|5.3|l results from the necessity of taking into account 
the rescattering diagrams in the axial channel of the pseudoscalar (w or rj') mesons (Fig. Q2Jl. The correlators Ij5.1|) are 
defined in such a way that the perturbative contributions are subtracted, 



n„ p (Q 2 ) = n tot (q 2 ) - n pert (q 2 



(5.5) 



The perturbative part is obtained from the non-perturbative part by simply setting the dynamical quark mass M (k) 
to zero. We extract the longitudinal and transverse parts of the correlators through the use of the projectors 



pL = t^_ 

l±v 2 : 



We first consider the transverse part of the V correlator, with the result 



K^(Q 2 ) = 2N C J 



S%(Q 2 ) = AN C 



d A k 



1 



(2^) 4 D+D- 

M + - M_ 



M+M. 



2 

k — k 1 

+ 3 x 



up 

M 2 - M 2 



,2 ,2 j (fe+fc--M + M_) 

n> 1 /v j -\- 



d4k M « V (fe) + l / J 



(2^) 4 !>(&) |/~ v "' ' 3fc2 _ ( fc + Q) 2 v 
where the notations k^ = k^ - Q> 1 {kQ)/Q 2 , k± = k± Q/2, 

M± = M{k±), D± = D(fc±) 



M'(fc) 



M (fc + Q) - M (k) 
(k + Q) 2 -k 2 



(5.6) 



(5.7) 



(5.8) 



(5.9) 





j/V'X/x/v, 4- 




J -1 /\s\s\ss* 4- *\ss/\s\j J -1 




FIG. 9: The resummation of the quark rescattering in the current-current correlator in the axial-vector channel. 




FIG. 10: The contact terms in the current- current correlators. 



have been introduced. The subtraction of the perturbative part amounts to the replacement 



1 



D+D- 



2 

/c_i_ k — k i 

+ 3 1 



np 



2 

k-i-k — k i 

+ 3 ± 



1 



1 



D+D- kik 2 



Further, we take the nonsinglet transverse projection of the A correlator and obtain 

d A k 1 



2N, 
3 k± 



(2*)' 









k-i-k — — k i 

. 3 . 



(* 



I -k 2 _ 



np 

^ (fc+fc _ + m + m_) - 



k\-k 2 _ 

2 



S$ (Q 2 ) = 5? (Q 2 ) - W c [ Al 

J (2ir 



d^k M(k)8 u2 (y/M(k + Q)-y/Mjk)) 



( 27r ) 4 D{k) 3 



(k + QY 
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(5.10) 



(5.11) 



(5.12) 



Let us consider the difference of the V and A correlators, where a number of cancellations takes place and the final 
result is quite simple, 



Q 2 Tl v - A (Q 2 ) 



AN, 



d i k 1 f 4 , 

j — 7~^T\ r^rr 1 M+M_ + -kj 

(2tt) 4 D (k 2 + ) D (k 2 _) \ 3 JL 



, M + - M_ 



A;^ - A; 2 



(5.13) 



'M-k. 



k 2 - k 2 



One may explicitly verify that the integrand of the above expression is indeed positive-definite, irrespectively of the 
choice of the mass function M(p). Thus the Witten inequality 12. 5|) is indeed fulfilled. 
At Q 2 = one gets the results consistent with the first Weinberg sum rule, 

- q 2 u v ' t (q 2 = o) = o, -g 2 n^ T (g 2 = o) = -/ 2 , -q 2 ii v ^ t (q 2 = o) = / 2 , (5.w) 

where the explicit definition of the pion decay constant (|4.26|) is used. This serves as a useful algebraic check. 

VI. MODEL PARAMETERS 



The parameters of the model are fixed in a way typical for effective low-energy quark models. We demand that the 
pion decay constant f n , (|4.26() . and the quark condensate (for a single flavor), < qq >, given by 



(6.1) 
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acquire their physical values. For simplicity, we take profile for the dynamical quark mass in a Gaussian form 

M(u) = M q exp (-2u/A 2 ) . (6.2) 

With the model parameters 

M q = 0.3 GeV, A = 1.085 GeV, (6.3) 

one obtains 

U = 93 MeV, (qq) = - (224 MeV) 3 , (6.4) 
where the quark condensate is supposed to be normalized at the scale of a few hundred MeV. 

VII. LARGE- Q 2 EXPANSION 

At large Q 2 one finds the following asymptotic expansion for the difference and sum of the correlation functions 
in the inverse powers of Q 2 (we do not display here the exponentially-suppressed terms coming from powers of the 
dynamical quark-mass form factor): 

■0(7^.(7.2) 



-Q 2 I% +A (Q 2 -> 00) = 2^r [ du- V 



4tt 2 J D(u) 



M{u)M'{u) 1- — --— ^ 



Q 2 3Q*J 6 Q 



The effective model considered here is designed to describe low energy physics. At high energies it is certainly not 
expected to reproduce all the details of the asymptotic standard operator product expansion of QCD. On other hand, 
it is possible that the OPE works well only at very short distances while the effective model is applicable at large and 
intermediate distances. With this hope in mind we proceed to analyzing the large-Q 2 expansions of the correlators, 
comparing them numerically to the OPE results, and trying to match the two approaches. It is important to note 
that the power corrections in the expansions i|7.1|l and (|7.2|1 have the same inverse powers of Q 2 as the OPE. 

We may now compare the expansion of the model correlators to the OPE , Eqs. Q2.6l2.7j l. In Eq. l|7.1j) the formally 
leading d — 4 term absent in the chiral limit in accordance with the second Weinberg sum rule and the OPE QCD. 
The second term in Eq. (|7.1|) (and the last term in Eq. (|7.2|) ') is proportional to the derivative of the gluon condensate, 
and via equations of motion it reduces to the four-quark condensate term appearing in the OPE, Eqs. i|2.6l2.7|) . Let 
us compare the numerical estimates for the local d — 6 terms obtained from the QCD sum rules and from the nonlocal 
chiral quark model, labeled as NxQM: 

[0%- A ] QCDsr « -(1 - 2) • 10- 3 GeV 6 , (7.3) 
[0%- A ]° PET = -(3.4 ±1.1)- IO- 3 GeV 6 , 
[0%- A f xQM = -1.1- HT 3 GeV 6 . 

The first estimate is found on the basis of low energy theorems and QCD sum rules 0, while the second estimate 
is made with the help of the r-decay data ^l| . The result of the present model is closer to the standard estimate 
obtained from the low-energy phenomenology. Similar features of the short-distance behavior of the correlators were 
found in the instanton model |13| . 

In the V + A correlator Ij7.2jl the short-distance expansion contains, in addition to the contributions coming from 
the local operators, the unconventional terms originating from the nonlocal operators of dimension d = 2, 4 and 6 
(the first terms in Eq. 17. 2(1 *1. This kind of unconventional terms has recently attracted attention due to the revision 
of the standard OPE as well as the lattice findings, where the unconventional power corrections in the vector 
correlators were reported j^] • The appearance of this correction is usually related to the existence of the lowest d = 2 

condensate V which is due to an apparent gauge non-invariance, absent in the standard OPE. However, in 

a series of papers (|H(], 001 , an d references therein) it was argued that it is possible to define the nonlocal operator 
with the lowest dimension in a gauge-invariant way. This situation is very similar to the famous spin-crisis problem 
(c/. H2). Analogously, in that case there is no twist- two gluonic operator that may contribute to the singlet axial 
current matrix element, yet, it is possible to construct the matrix element from nonlocal operators [Bljj ]. We thus 
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see that our effective nonlocal model shares these unusual effects generated by the internal nonlocalities of the quark 
interaction. Furthermore, the lowest-dimension power corrections are naturally present in the approaches similar to 
the analytic perturbative theory j^. In that case in order to compensate the effects of the ghost pole in the strong 
coupling constant, the d = 2 power term is added. As discussed in Rcf. ji^, the justification of the appearance of the 
unconventional power corrections at the same time means that the standard OPE is valid only at very large momenta. 

We also wish to comment that in the model expansion of the V + A correlator there are no explicit terms with 
the gluon condensate of dimension d — 4. The appearance of this term in the nonlocal chiral quark model would 
corresponded to the local matrix element 

iV c f du uM 2 (u) 



At: 2 J D(u) ' 

that is related to the gluon condensate through the gap equation 1)3. 7f) [44|. However, the coefficient of this term is 
equal to zero. This is due to the simple form of the quark propagator (|3.6I) . which docs not allow gluonic correlations 
between different quark lines. The similar situation occurs in the QCD sum rules calculations (in the fixed point 
gauge), where nonzero contribution comes from the diagram with quark lines correlated by soft gluon exchange. 
These (numerically small) correlation terms may be reconstructed in the effective model by introducing the gluonic 
field in the effective action (|3.2|) by gauging kinetic and interaction terms (see also cf. (34]). From other hand the 
d = 4 term appears in (|7.2(l as a nonlocal matrix element. 

We end the discussion of the short-range behavior of the correlators by giving the numerical estimates of the 
additional terms appearing in Eq. I|7.2[) : 

[0X +A t:i: = -£/ du^M(u)M' { u) = 5.0- 10- GeV 2 , (7.4) 

[OVXI: = ^/d«^M(„)M' („) = -1.8-10-3 GeV 4 , 

r V+A^XQM = ANc f du ^_ M{u)M ' {u) = _ T Q. l0 -4 GeY 6 
L 6 nonloc 3 4^2 J £) W V I 

The sum of these terms, taken in the interval of momenta Q 2 ~ (1 -j- 2) GeV 2 where the model large-Q 2 expansion is 
expected to be valid, agrees reasonably well with the coefficient of the leading power correction in Eq. l|2.7[l 

[0% +A ] QCDst = 3.0-lQ- 3 GeV 2 , (7.5) 

where we have taken the estimate (a s /7r) A 2 = —0.12 GeV 2 from Ref. |42^ . 

Through the use of the factorization hypothesis l|2.9|) it is predicted that the chirality flip matrix element 0\~ A is 
strongly enhanced in absolute value over the chirality nonflip one Oq + 

[0%- A /0% +A ] lactoI = -4.5. (7.6) 

In the nonlocal chiral quark model the chirality flip matrix element 0§~ A is given by the local matrix element, but 
the chirality nonflip one is a mixture of the local and nonlocal matrix elements which transform to each other 

by integration by parts. We find that their ratio 

[Oj- A /OX +A f xQM «-3.2 (7.7) 

has the same tendency as predicted in (|7.6(l . It happens due to partial compensation of contributions of the local and 
nonlocal matrix elements into Oq +a . 



VIII. LOW-ENERGY OBSERVABLES AND THE ALEPH DATA 

Let us now consider the low-energy region where the effective model (|3.2|) should be fully predictive. From (|5.13|) 
and the DGMLY sum rule (|2.4(l we estimate the electromagnetic pion mass difference to be 

[m n ± - m n o] NxQM = 4.2 MeV, (8.1) 
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which is in remarkable agreement with the experimental value (after subtracting the — m u effect) [(j 

K± -"V>] exp = 4.43 ±0.03 MeV. (8.2) 



It is interesting to estimate the electric polarizability of the charged pions, |57l l58j . With the help of the DMO sum 
rule |73| . as done by Gerasimov in |56| . we find 



a 



3 

where Idmo is the left-hand side of the DMO sum rule l|2.3j) 



a 7T± = 



L DMO 



1 f s ° ds 

Idmo (s ) = ~ N ( s ) _ a i ( s )l • 



(8.3) 



(8.4) 



Equation l|8.3|l can be interpreted as a sum of the center-of-mass recoil contribution and the intrinsic pion polariz- 
ability. In 56] it was demonstrated that in model calculations there occurs a delicate cancellation between the two 
contributions of Eq. I|8.3(l . This requires the calculation of both terms consistently at the same level of approximations. 

With the experimental value for the pion mean squared radius |3S| and the value of the Idmo integral estimated 
from the ALEPH and OPAL data 



[Idmo K)] cxp = (26.3 ± 1.8) • 10" 3 (8.5) 



one gets from Eq. (|8.3|) the result Q 



[af±] cxp = (2.71 ± 0.88) • lO" 4 fm 3 . (8.6) 



From Eq. (|2.3|l also follows the relation obtained by Terentyev [^J, which relates the pion polarizability and the 
pion axial- vector form factor, 



<& = ^% (8-7) 



The last relation, used with the known values for Fa, yields 

a%± = (2.69 ± 0.37) • 10~ 4 fm 3 , (8.8) 

which is very close to (|8.6II . 

Let us estimate the electric polarizability of the charged pions within the nonlocal chiral quark model. By calculating 
the derivative of II^ - " 4 (Q 2 ) at zero momentum we estimate the left hand side of the DMO sum rule as 

[Idmo (*o - ™)}n x qm = ^.2 ■ 10" 3 . (8.9) 
The value of the pion charge radius squared, 

[<^>Wm= - 33 fm2 > ( 8 - 10 ) 

obtained in our model from the derivative of the charged pion form factor, is close to its limit of the local chiral model, 
found by Gerasimov long ago |H^. 

[<^>W=4^2 = - 34 fm2 - ( 8 - n ) 

The model predictions for (r 2 .) and Idmo are somewhat smaller than the experimental values given above. The 
reason for this discrepancy may be attributed to vector meson degrees of freedom, neglected in our treatment, and 
to pion loops absent in the large-JV c limit. However, these unconsidered contributions are essentially canceled in the 
combination (|1.4|l defining the electric pion polarizability (see, e.g. 0] for discussions). From H8.3(l we find with 
values given in Eqs. fSJJl and (|CT7|) the value 

K4 x QM = 2 ' 9 - irV . ( 8 - 12 ) 
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FIG. 11: Inclusive vector minus axial vector spectral function, v\ — ai, measured by the ALEPH collaboration 



which is close to experimental number (|8.6|) and also to the prediction of the chiral perturbation theory at the one-loop 
level El, 



X PT 



2.7 • HT'W 



.13) 



Let us note also that (|8.12|) is a factor of 2 smaller from the estimates obtained in a local chiral quark model [6l| . 
a^ ± — 5.8 ■ 10~ 4 fm 3 . We thus see that the model prediction for the pion polarizability, Eq. I|8.12|l . is in a very 
reasonable agreement with the experimental data. 

Next, we compare the model correlators with the ALEPH data, presented in Fig. The ALEPH and OPAL data 
integrated up to the r mass satisfy all chiral sum rules within the exp erimental uncertainty, but the central values 



expe 

differ significantly from the chiral model predictions. Following Ref. we use Sq = 2.5 GeV 2 as an upper integration 
limit, the value at which all chiral sum rules are satisfied assuming that Ui(s) — a>i(s) = above sq. Finally, a kinematic 
pole at q 2 = is added to the axial-vector spectral function. The resulting unsubtracted dispersion relation between 
the measured spectral densities and the correlation functions becomes 



1 



ax(s) ft 



Q 2 



where f 2 is given by the WSR I, 



1 f s ° 



(8.14) 



(8.15) 



Having transformed the data into the Euclidean space, we may now proceed with the comparison to the model, which 
obviously applies to the Euclidean domain only. Admittedly, in the Euclidean presentation of the data the detailed 
resonance structure corresponding to the p and a\ mesons seen in the Minkowski region is smoothed out, hence the 
verification of model results is not as stringent as would be directly in the Minkowski space. In Fig. 1121 we show the 
normalized curves corresponding to the experimental data and the model prediction. We also show the prediction of 
the model of Ref. ^5| (minimal hadronic approximation, MHA [62l l63lp 
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FIG. 12: Euclidean-momentum correlation function, — Q 2 TIv-a(Q 2 )/ f%, constructed in the present model (solid line), in the 
model of Ref. [l5l| (dashed line), and reconstructed via Eq. 18.141 from the ALEPH experimental spectral function of Fig. 1111 
(dash-dotted line). 

where the contributions of the p and a\ mesons are taken into account with the model parameters M p — 0.750 GeV 
and f a = M 2 /M 2 = 0.5. Other parameters are constrained by the Weinberg sum rules. As demonstrated in Ref. [lRj . 
the good agreement between the data and the model predictions is far from trivial, since many analytic approaches 
discussed in the literature meet definite difficulties in the description of the ALEPH data in the region of moderately 
large Q 2 . In Fig. (|13|) we also present the ratio of the nonperturbative parts of the V — A l|5.13[) to V + A correlators 
in the nonsinglet channel. 

To conclude this Section we wish to recall that quite similar calculations of the vector and isovector axial-vector 
correlators within a nonlocal model were done some time ago by Holdom and Lewis . There are certain differences 
in the form of the nonlocal interaction and, as a consequence, the form of quark-current vertices is different. A more 
principal difference is that the authors of |4fij have chosen a "two phase" strategy, describing the low-energy part 
of correlators by full nonperturbative vertices and propagators, while the high energy parts were computed in the 
approximation when one of the vertices is local. In this case the problem of matching of two regimes occurs already 
at rather low energy scales. In the present calculations one prolongs the applicability of the model up to moderately 
large energies, which inter alias results in a good description of the ALEPH data. On other side, we have to admit 
that one of the goals of both approaches, namely the finding of a direct correspondence between the effective model 
calculations and the OPE in QCD has not yet been reached (see Sect. VII). To make correspondence closer it is 
necessary to supply the model with a more detailed information on the soft quark-gluon interaction |74|. 



In this Section we demonstrate explicitly the transverse character of the V and isovector (IV) A correlators (Figs. 
||}5J and I|1U|)1. For the longitudinal component of the V correlator we get 



IX. CURRENT-CURRENT CORRELATORS (LONGITUDINAL PARTS) 




(9.1) 



and therefore 




(9.2) 



as it certainly should be by the requirement of the vector current conservation. 
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FIG. 13: The ratio of the nonperturbative parts of the V — A to V + A correlators in the nonsinglet channel. 



Further, we consider the longitudinal projection of the A correlator. Then, we get contributions from the one- 
quark-loop diagram 



Ki,i (Q 2 

the two-loop diagram in the isovector channel 

8N, 



AN C f d A k M(k) 



Q 2 J (2n) 4 D(k) 



[M (k + Q) + M (k)} 



the one-loop contact diagram 



d k M (k) nrrr, ^tttttt 

QtfDU VM{k + Q)Mik) 



-i 2 r 



d 4 k M 2 (k) 



-i -l 



(2irf D(k) 



AN C f d 4 k M(k) 



Q 2 J (2tt) 4 D (k) 



3M (k) +M(k + Q)- Ay/M (k) M (k + Q) 



and, finally, from the two-loop contact diagram in the isovector channel 

d A k M{k) 



st,2 (Q 2 



8N C 



(2tt) 4 D(k) 



M (k) - ^M(k + Q)M(k) 



d A k M 2 (k) 
(2tt) 4 D{k) 



The sum of all these contributions leads to the desired result 



Q 2 nf v (o 2 ) = o 



(9.3) 



(9.4) 



(9.5) 



(9.6) 



(9.7) 



consistent with the isovector axial current conservation in the strict chiral limit. 



X. SINGLET AXIAL VECTOR CURRENT CORRELATOR AND THE TOPOLOGICAL 

SUSCEPTIBILITY 

The cancellations in the longitudinal channels are consequences of the current conservation and follow simply from 
the application of the nonanomalous Ward-Takahashi identities. We have explicitly demonstrated this in the previous 
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Section in order to show the consistency of our calculations. The issue becomes important when we consider the 
longitudinal part of the isosinglet axial- vector current correlator which is not conserved due to the Ua (1) axial Adlcr- 
Bell-Jackiw anomaly. This channel is dominated not by the pion, but by the //-meson intermediate state. Thus, 
in addition to the nonlocal interaction present in Eq. 13.2(1 we also need to include the interaction 1(3.5(1 . where an 
exchange of the "77" singlet meson, the SU (2) analog of the 7/ meson, occurs. 

It is well known that due to the anomaly the singlet axial- vector current is not conserved even in the chiral limit, 

P J*° (x) = 2N f Q 5 (x) , (10.1) 

where Q 5 (x) is the topological charge density. In QCD it is defined as Q 5 (x) = (a a /8n)G*(x)G"(x), where is 

the gluonic field strength, and G a is its dual, G^„ = e^vXaG^. The correlator of the singlet axial- vector currents 
has the same Lorentz structure as in Eq. 1(1.2(1 . but without flavor indices and with T a = 1. In the chiral limit its 
longitudinal part is related to the topological susceptibility, the correlator of the topological charge densities Q5 (x) , 



X{Q 2 )=^ / d 4 xe^(0|T{Q 5 (x)Q 5 (0)}|0), (10.2) 

by the relation (see, e.g., |b3 |1 



(Q 2 ) = Wt x (Q 2 ) . (10.3) 



At high Q 2 the OPE for X (Q 2 ) predicts \(£ 



X (Q 2 - 00) = (G^) 2 ) + 0(Q- 2 ) + 0(e-n (10.4) 

where the perturbative contributions have been subtracted, and the exponential corrections are due to nonlocal 
instanton interactions |64| . 

At low Q 2 the topological susceptibility \ (Q 2 ) can t> e represented as a sum of contributions coming purely from 
QCD and from hadronic resonances, (6^|. Crewther proved the theorem [(Jl) that x (0) =0 in any theory where at 
least one massless quark exists (the dependence of x (0) on current quark masses has been found in |S3|Z3). Also, the 
contributions of nonsinglet hadron resonances are absent in the chiral limit. Thus, in the low-Q 2 limit for massless 
current quarks one has 

[* O? 2 - °)] xQ cd = -QV(0) + 0(Q 4 ). (10.5) 

The estimates of x'(0) existing in the literature are rather controversial |7q : 

x'(0) = (48 ±6 MeV) 2 [69], x'(0) = (26 ± 4 MeV) 2 [70]. (10.6) 

This makes further model estimates valuable. 

Now we turn to the model calculations. The bare isosinglet axial- vector current obtained from the interaction terms 
((3.3(1 and ((3.5(1 by the rules described in Sect. IV becomes 

rf(k,q,k' = k + q) = w - 75 (k + k')^ ^— p J -+ (10.7) 



75%2VM(fc')M(fc) 



<1 2 



where JAp(q 2 )is defined in Eq. ((4.7(1 . In order to get the full current we have to consider rescattering in the channel 
with the quantum numbers of the singlet pseudoscalar meson, "77" , which results in 

(\/M(F)- ^M(k)Y 

Tf(k, q,k'=k + q) = 7m75 - 7 5 (fc + fc %" ^— p '— (10.8) 
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Because of the singlet axial anomaly this current does not contain the massless pole anymore, since according to Eq. 
(|4.27|) one has at zero momentum: 



l-GJpp(q 2 



G 



A 

Mf 



(10.9) 



Instead, it develops a pole at the "77" meson mass. So, within the model considered the same mechanism is responsible 
for the sponteneous breaking of chiral symmetry and violation of the ?7a(1) symmetry. 
The vertices satisfy the anomalous Ward-Takahashi identities: 



q„Tf(k, q, k' = k + q)= q l5 - 75 W(k') + M(k)} + l5 2y/ M(k')M(k) ( 1 - — + G 



a 



„JAp(q 2 ) 



Ml 



and 



q^f(k, q,k' = k + q) = q l5 - 75 [M(fc') + M(k)} + 75 



2y/M(k')M(k) 
l-G'J PP {q 2 ) 



1 - 



G' 



G 



(10.10) 



(10.11) 



where the last term is due to the anomaly. Thus the QCD pseudoscalar gluonium operator is interpolated by the 
pseudoscalar effective quark field operator with coefficient expressed in terms of dynamical quark mass. In the 
effective quark model the connection between quark and integrated gluon degrees of freedom is fixed by the gap 
equation (|3.7|l 44]. By considering the forward matrix element (q = 0) one deduces that the singlet axial constant 

is not renormalized within our scheme: (Q) = 1. In order to get reduction of the singlet axial constant ("spin 
crisis") we need to consider the effects of polarization of topologically neutral vacuum (see, e.g., 

It is instructive first to consider the longitudinal part of the correlator of the local vertex, 7^75, and the nonlocal 
vertex of Eq. (|10.8|) , which is the construction of Pagels and Stokar 48] . In this model the decay constant is defined 
by 



U.PS - ^2 



du 



2M(u) 2 - uM(u)M'(u) 
1 Diju) • 



(10.12) 



Then, through the use of Eq. (|10.3(l we get for the topological susceptibility the result 

G-G' 



(2N f ) 2 X PS (Q 2 ) 



N c Nj 

G[l-G'J P p(Q 2 )} 4tt 2 
with the coefficients of the low-Q 2 expansion given by 



Xps (0) = 0, x'ps (0) = 



2N 



f 



(10.13) 



(10.14) 



Hence, the result is consistent with the Crewther theorem and it provides the estimate of x'ps (0) ~ (39 MeV) 
obtained for Nt = 3. The second relation may be rewritten in the form resembling the generalized Goldberger- 
Treiman relation, as advocated by Veneziano and Shore |7Cj . Indeed, by using the standard Goldberger-Treiman 
relation, (|4.24|) . valid in a given model, one finds 



M q =gZ S yf2N fX ' PS (0), 



(10.15) 



which is just the quark-level relation from Ref. [jjj- At large Q 2 the quantity xps (Q 2 ) decreases according to the 
power of the dynamical quark form factor. 

Now let us turn to the full model calculation. Proceeding in a similar manner as in the previous Sections we get 
the topological susceptibility in the form 



(2N f f X (Q 2 



G' 

2N f [l~ — 



M 2 J PP (Q 2 ) 1 



Q 2 JttA (Q 2 
G' 



G'J A p(Q 2 ) 



1 



Mf 



Ml 



Jap(Q 2 ) 



GJap(Q 2 
M 2 



G 

Ml 



4N c N f 



d 4 k M(k) 
(2tt) 4 D(k) 



l-G'J P p(Q 2 ) 

G-G' 
G[1-G'J PP {Q 2 
2 " 



(10.16) 



M (k) - y/M (k + Q)M(k) 
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where we have used the relation between the singlet current correlator and the topological susceptibility of Eq. i|10.3fl • 
At large Q 2 one obtains the power-like behavior consistent with the OPE prediction l|lU.4fl . namely 

9,0 s 2N f M 2 ( G'\ 
- (27V/) 2 X {Q 2 - 00) = —^L (l-^y (10-17) 

At zero momentum the topological susceptibility is zero 

X (0) = 0, (10.18) 
in accordance with the Crewther theorem. For the first moment of the topological susceptibility we obtain 

x'(0)^|/ 2 (2-§) + (l-§) 2 ^(0)}, (10.19) 

where f 2 and J' AP (0) are defined in Eqs. (|4.2ti|) and (|4.9|) . respectively. 

In order to get numerical results we need to specify further the details of the model. We consider two possibilities. 
One involves the interaction with the exact symmetry as provided by the 't Hooft determinant, G = — G. For the 
second possibility the constants G and G' are considered as independent of each other, and their values are fixed with 
the help of the meson spectrum. In this more realistic scenario one has approximately the relation G' = 0.1 G (for 
typical sets of parameters, c.f. Ref. |23j). Then we get the following estimates for the first moment of the topological 
susceptibility: 

x'(0) = (55 MeV) 2 {G 1 = -G), (10.20) 
x'(0) = (50 MeV) 2 (G' = 0.1 G). (10.21) 

For the above estimates we have taken Nf — 3. Since the flavor number enters only through the factor of 2Nf present 
in the definition (|10.3|) . in this sense it is external to the model and its inclusion is very simple. We can see that the 
model gives the values of x'(0) which are close to the estimate of Ref. |29j. In Fig. 3 we present the model predictions 
for the topological susceptibility at low and moderate values of Q 2 for the cases of the full l|10.16fl and Pagels-Stokar 
(|10.13fl model calculations. 

We should note that the predictions of our model have a limited character because we have used the SU / (2) model 
in the chiral limit and have not considered mixing effects. However, our final result is formulated in terms of a 
physical observable, f T , and thus we believe that the presented predictions may be not far from more realistic model 
calculations. In the region of small and intermediate momenta our results are quantitatively close to the predictions 
of the QCD sum rules with the instanton effects included (6^|. In our opinion, both the instanton-based calculations 
(our model with (G 1 = — G) and the interpolation of the model |U) overestimate the instanton contributions in the 
region Q 2 ~ 0.5 — 2 GeV 2 . It would be interesting to verify the predictions given in Fig. El by the modern lattice 
simulations. 



XI. CONCLUSIONS 



In this work we have analyzed the nonperturbative parts of the Euclidean-momentum correlation functions of the 
vector and axial-vector currents within an effective nonlocal chiral quark model. To this end, we have derived the 
conserved vector and isotriplet axial- vector currents and demonstrated explicitly the absence of longitudinal parts in 
the V and nonsinglet A correlators, which is consequence of the gauge invariance of the present approach. On the 
other hand, the singlet A correlator gains an anomalous contribution. From the properties of the V — A correlator we 
have shown the fulfillment of the low-energy relations. The values of the 7r ± — tt° electromagnetic mass difference and 
the electric pion polarizability are estimated and found to be remarkably close to the experimental values. In the high- 
energy region the relation to OPE has been discussed. In particular, the estimate of the 1/Q 2 coefficient is in agreement 
with the recent lattice findings and the modified OPE phenomenology. We stress that the momentum dependence 
of the dynamical quark mass is crucial for the fulfillment of the second Weinberg sum rule. The combination V — A 
receives no contribution from perturbative effects and provides a clean probe for chiral symmetry breaking and a test 
ground for model verification. We have found that our model describes well the transformed data of the ALEPH 
collaboration on the hadronic r decay. The combination V + A, on the other hand, is dominated by perturbative 
contributions which are subtracted from our analysis. By considering the correlator of the singlet axial- vector currents 
the topological susceptibility has been found as a function of the momentum, and its first moment is predicted. In 
addition, the fulfillment of the Crewther theorem has been demonstrated. 
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FIG. 14: Topological susceptibility versus Q 2 predicted by the model with G' — 0.1 G, Eq. 110.161 . (solid line), and by the 
Pagels-Stokar construction, Eq. I1U.13H , (dashed line). 
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[73] In xPT in the one-loop approximation the right-hand-side of the DMO sum rule is expressed through the low-energy 

constant Lio- The extraction of this constant from the experiment was considered in 
[74] Let us also refer to other important works |46|. |47|. where the problem of connection between the effective 4-fermion 

models and QCD has been discussed. 
[75] Consistencies of the axial- vector current conservation and Ua(X) problem is also discussed in [oSf. 
[76] The results obtained in Ref. |7l| concerning x(0) and x'(0) contradict to the low-energy theorems. 



